Phase separation in optical lattices in a spin-dependent external potential 
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We investigate the phase separation in one-dimensional Fermi gases on optical lattices. The den- 
sity distributions and the magnetization are calculated by means of density-matrix renormalization 
method. The phase separation between spin-up and spin-down atoms is induced by the interplay of 
the spin-dependent harmonic confinement and the strong repulsive interaction between intercom- 
ponent fermions. We find the existence of a critical repulsive interaction strength above which the 
phase separation evolves. By increasing the trap imbalance, the composite phase of Mott-insulating 
core is changed into the one of ferromagnetic insulating core, which is incompressible and originates 
from the Pauli exclusion principle. 
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I. INTRODUCTION 

Ultracold atoms in optical lattices provide a new test 
bed for interacting quantum many-body systems 
Fermionic atoms in optical lattices can be used to re- 
alize the clean Fermi-Hubbard model, which is free of 
lattice defects, impurities, and phonons, in contrast to 
those in solid-state systems. Over the past few years, 
many interesting phenomena were observed in optical lat- 
tices, for example, the Fermi surface of the atoms in the 
lattice, the transform from a normal state into a band 
insulator 0, and fermionic superfluidity of attractively 
interacting fermions Two other major breakthroughs 
achieved recently in fermionic superfluid are: the BEC- 
BCS crossover [4] and imbalanced superfluidity 

The spatial inhomogeneity due to the confinement 
essential for ultracold atomic experiments is always 
present, which leads to a spatially varying local den- 
sity distribution and normally invalidates a reliable an- 
alytical method usually used in the homogeneous sys- 
tem. Many numerical schemes such_a_s_the density-matrix 

uantum Monte 
l3l , and density- 



renormalization group (DMRG) 
Carlo [9I4T2I , exact diagnolization 
functional theory based on the exact Bethe-ansatz so- 
lution ^] are used in studying the many interesting 
quantum effects in spin-balanced or imbalanced systems. 
Among them, intriguing properties such as phase sepa- 
ration in a trap and the transition from superfluidity to 
a normal state have attracted a great deal of attention 
both experimentally and theoretically. 

In the experiments, a phase separation was observed 
between the normal component and the superfluidity of 
interacting fermionic atom gases with imbalanced spin 
populations Q. In theory, the mean-fleld approach pro- 
vides a qualitative explanation of the phase separation 
of imbalanced fermionic atom gases in a trap 15]. The 
imbalance of the two species with A^-|- ^ can be 




produced by different trapping frequencies [16|, namely, 
spin-dependent trapping potentials. Phase separation 
can occur in trapped spinor boson gases with a weak 
anisotropic spin-spin interaction [l7| and in multicom- 
ponent Fermi gases with different values of the scatter- 
ing lengths and particle number [l^ . In two-dimensional 
optical lattices, the phase separation due to the imbal- 
anced mixture, antiferromagnetic order [Tq! ]. and pairing 
symmetries [2Q!] is investigated by the mean-field theory. 
For a one-dimensional (ID) system of two-component 
Fermi gases in a continuous space, it is found that there 
exists a critical interaction strength beyond which one 
atomic component expels another from the center of the 
trap [21]. For a ID attractive Hubbard model, a phase 
separation between the condensate and unpaired major- 
ity atoms occurs for a certain range of the interaction 
and population imbalance. At T = beyond a critical 
spin polarization, the phase separation always exists no 
matter how strong the interaction is j22]. For a ID re- 
pulsive Hubbard model, the phase separation due to the 
different trap frequencies is discussed within the local 
magnetization by the spin-dependent density-functional 
theory @]. 

In the present work, we are interested in the phase 
separation between different fermion species induced by 
the spin-dependent external potentials. The interplay 
between the external spin-dependent potentials and the 
repulsive interaction of intercomponent fermions will be 
explored. 



II. THE MODEL 
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We consider a two-component Fermi gas in a tube with 
Nf atoms and Ns lattice sites with the unit lattice con- 
stant, which can be described by a one-band inhomoge- 
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neous Fermi-Hubbard model [24 
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where the spin degrees of freedom a =t, ]. are pseudospin- 
1/2 labels for two internal hyperfine states and Ci^ (c|g.) 
are fermionic operators annihilating (creating) particles 
with spin a in a Wannier state at site i. fii — fiiu — 
So- ^LciCT is the total site occupation operator, t is 
the tunneling between the nearest neighbors, U is the 
strength of the on-site interaction, and Va describes the 
strength of the spin-dependent harmonic trapping poten- 

tials Fhar,. ^V^[i- {Ns - l)/2]'. 

The inhomogeneous Fermi-Hubbard model can be real- 
ized by a strong confinement in transverse directions [25j 
with an additional periodic potential applied along the 
tube. Concerning the experimental realization of the 
spin-dependent external potentials, one can use magnet- 
ically trapped Fermi mixtures of a. particular atom in 
the two different hyperfine states [26|, l27ll. or two dif- 
ferent trapped atoms of unequal masses [28|, where the 
different magnetic moments make V^. In the ex- 

periment of two ^"^K fermion species, the ratio of fre- 
quencies Vf/Vi = \/9/7 is discussed [2^. In optical 
lattices, a spin-dependent optical trap can be realized 
by asymmetrically detuning the laser frequencies with 
respect to the two hyperfine states [l^l- Experimen- 
tally, the atomic density we calculated is the most conve- 
nient and clear observable detectable by electron beams, 
high-resolution cameras, or noise interference. Recently, 
a composite phase of an incompressible Mott-insulator 
phase in the core was identified [30| . where the core is 
composed of strongly repulsive fermionic atoms in two 
hyperfine states. It is shown how the system evolves by 
increasing confinement from a compressible dilute metal 
into a band-insulating state, which also provides a way to 
polarize a spin-balanced system where — — Nf/2 

The homogeneous ID Fermi-Hubbard model belongs 
to the universality class of Luttinger liquids. At zero 
temperature, the properties of this model in the thermo- 
dynamic limit {Na,Ns oo, but with finite Ncr/Ng) are 
determined by the fillings no- — N^/Ns and by the dimen- 
sionless coupling constant u = U/t. According to Lieb 
and Wu ^T], the ground state (GS) properties for dif- 
ferent fillings in the thermodynamic limit are described 
by the coupled integral equations (for details see Refs. 
[32,33]). 

For the inhomogeneous system described by Eq. ([l}, 
the coexisted phases induced by the external spin- 
independent trapping potentials (Vj_ = Vj,) were well 
identified by many authors 0, [13, Hj-Is^ . We focus in 
this work on the spin-dependent potentials (V^ 7^ Vj,) by 
applying the DMRG techniques, performed by using the 
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FIG. 1: (Color online) Schematic illustrations for spin- 
dependent harmonic potentials Viiar.o- (in units of t) and the 
density distributions rio- (in units of the lattice constant) of 
both spin-up and spin-down atoms in the presence of inter- 
actions. The left panel (a) is for the system of small trap 
imbalance, where the spin-up and spin-down atom mixture 
in the center forms a PM region and the right panel (b) for 
the system of large trap imbalance, where the PS region is 
formed. 



ALPS libraries [37|. During our DMRG calculations, the 
states kept are 500 to 1000 so that we can restrict the 
cut error to be less than 10~^^. 



III. NUMERICAL RESULTS 

In this section we present our numerical results. In the 
following discussion we keep the total number of parti- 
cles constant f = 40) and vary the number of spin-up 
and spin-down atoms in the system. We characterize the 
confinement imbalance by defining the ratio between the 
spin-up and spin-down dependent external potentials as 



(2) 



In Fig. [T] we show schematic plots for the spin- 
dependent harmonic potentials and the density distribu- 
tions of the two fermion species with small or large con- 
finement imbalances. The effects of 7 are manifested in 
that two atoms coexist for small 7 where the spin-up and 
spin-down atom mixture in the center of the trap forms 
phase mixing (PM) region and separate with only spin- 
up atoms left for large 7 where the phase separated (PS) 
region is formed. The PS region is determined with the 
local occupation in the trap center (i.e., i=0) satisfying 
«o < 10"'^. We distinguish in the following the different 
phases by showing the atomic density profiles and the lo- 
cal magnetization for different repulsive interactions and 
confining strengths in the system of spin-unpolarized or 
spin-polarized atoms. 

First, we study the phase separation between the 
two-component fermions induced by the interplay be- 
tween the repulsive interaction and the spin-dependent 
parabolic potentials in the unpolarized system of an equal 
number of spin-up and spin-down atoms (N^^ — — 20) 
and Vj, = 1.0 x 10"'^. The lattice size chosen here and 
in the following is always large enough to make sure that 
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FIG. 2: (Color online) Phase diagram as a function of u and 
the confinement ratio 7 determined by the DMRG technique. 
The system consists of N-f = 20 and A*'^ = 20 fermions. The 
spin-down trap strength is V4, = 1.0 x 10~^. The arrow in 
the top indicates the position where 7=1 and the arrow in 
the right where Uc — 1.64. The solid line is a power-law fit 
u = Uc + a/{'y — 1) to the data with a = 10.932. The two 
phases PM and PS are manifested in Fig. [T]and explained in 
the text. 



the GS densities smoothly drop down at the edges. In 
Fig. [2I the phase diagram is shown as a function of u 
and the confinement ratio 7. Two regions are seen: the 
PM region with both spin-up and spin-down mixtures in 
the center of the trap and the PS region with only spin- 
up atoms remaining in the center. A critical interaction 
strength Uc = 1-64 is obtained, below which there is no 
phase separation no matter how large the confinement 
ratio. For the system considered here, the condition in 
which the phase separation happens can be simply fit- 
ted by a power-law relation u = Uc + a/{'y — 1), with 
a = 10.932. We further illustrate in Fig. [3] an explicit 
example by choosing 7 = 3 and changing the interac- 
tion strength. We confirm that there exists a critical 
value of the interaction strength {u = 8) beyond which 
the spin-down atoms are depleted from the center of the 
trap and repelled into the periphery regions between Inl- 
and V^. In this case, a phase separation begins to ap- 
pear (i.e., the Fermi components tend to stay in different 
spatial regions). Thus, upon approaching the phase sep- 
aration point and beyond, the local polarization of the 
atomic gases in the center becomes stronger and stronger. 
When the complete phase separation is realized, fermions 
become fully polarized due to the strong repulsive inter- 
action. As a result, spin-up atoms locate in the center 
and spin-down atoms at the periphery of the trap, which 
is clearly seen in Fig. ^c) for u = 20. Upon reaching 
the complete phase separation, further increasing the re- 
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FIG. 3: (Color online) Ground state density profiles for n;^, 
riii and local magnetization nn = (ni-|-— ni^)/2 as a function of 
i in the spin-dependent external potentials of 7 = 3. The spin- 
down trap strength is 10""^. The system consists of A'^ = 20 
and A^^ = 20 fermions. Three different interaction strengths 
are shown: (a) u = 1, (b) it = 8, and (c) u = 20. The local 
magnetization rrii is shown in (d). The solid line connecting 
the symbols serves as a guide for the eyes. For comparison, 
the GS densities of the noninteracting case (u = 0) for spin- 
up (bold solid line) and spin-down (thin solid line) atoms are 
included in (a), and the corresponding local magnetization 
(thin solid line) is also shown in (d). We find that the re- 
pulsive interaction can induce a complete phase separation 
between the two components in the spin-dependent external 
potentials. 



pulsive interaction only makes the spin-up density a little 
more confined and spin-down density more spread out. In 
Fig. El^d), we plot the local magnetization of the system, 
which is defined as rrii — (rti-f- — ni^)/2. For the strong 
repulsive interaction where the phase separation begins 
to evolve m.i changes from negative to positive with a big 
slope. 

Now, let us concentrate on the phase separation in- 
duced by spin-dependent parabolic potentials. In Figs. 
|4] and [5] we study the polarized systems of an unequal 
number of spin-up and spin-down particles (iV-|- = 30 
and Ni — 10) with weak {u = 1) and strong {u = 4) 
repulsive interactions. We illustrate the effects of the 
confinement ratio on the local density distributions and 
the local magnetization. We increase 7 by keeping the 
spin-down external potential as invariant and increas- 
ing V-f (i.e., 7 > 1). In Fig. SI the density profiles for the 
weak repulsive interaction {u — 1) are shown with dif- 
ferent confinement imbalances (7 = 1,3, and 6). While 
increasing the confinement for the spin-up atoms, the in- 
teraction between the spin-up and spin-down atoms in 
the center of the trap repels the spin-down atoms into 
the edges of the trap. However the repulsive interaction 
is not strong enough and only a small amount of phase 
separation appears. From Fig. [Sfb), we can see that, in 
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FIG. 4: (Color online) Ground state density profiles for 
n-i^, riii and local magnetization rrii as a function of i for 
the system of weak repulsive interaction (u = 1) in the 
spin-dependent external potentials. The system consists of 
= 30 and = 10 fermions with V; = 2.5 x 10"". Three 
different ratios of confining potentials are shown, (a) j — 1, 
(b) 7 = 3, and (c) 7 = 6. The local magnetization rrii is 
plotted in (d). From (d), we can see that upon reaching the 
phase separation point and beyond, the local magnetization 
rrii becomes more negative at the periphery and more posi- 
tive in the bulk region of the trap signaling that more spin- 
down fermions are repelled from the center and more spin-up 
fermions are constrained there. 
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FIG. 5; (Color online) Same as Fig. |4]but for the system of 
strong repulsive interaction (u = 4). From (b), we notice that 
a complete phase separation occurs. For comparison, in (a)- 
(c), the CS densities of the noninteracting case {u = 0) for the 
spin-up (bold solid line) and spin-down (thin solid line) atoms 
are also plotted. In (d) we include the corresponding local 
magnetizations for the noninteracting case with 7=1 (bold 
solid line), 3 (thin solid line), and 6 (dotted line), respectively. 



the system of strong repulsive interaction {u = 4) and a 
large trap imbalance (7 — 3), almost all the spin-down 
atoms are repelled from the bulk of the trap and a com- 
plete phase separation is realized. Due to the depletion of 
the spin-down fermions, fully polarized gases of spin-up 
fermions are obtained, as can be seen in Figs. EJb) and 
[5jc). For comparison, the GS density distributions of the 
spin-up and spin-down atoms for the noninteracting case 
(u — 0) are also included, where no phase separation is 
observed. We conclude that the intercomponent interac- 
tion is essential in achieving a phase separation between 
the two-component fermions in a spin-dependent trap. 

The local magnetization nii, in Figs. lUJd) andlS^d), 
gives another signature of the phase separation. For small 
7, a flat region of rrii is seen in the center of the trap and 
two bumps are shown at the edges with the excess spin- 
up atoms. The increase of the trap imbalance and the 
repulsive interaction strength shows a signature that rrii 
is more negative at the edges, that is, more and more 
spin-down atoms are repelled from the center of the trap 
and accumulate at the periphery region between Vf and 

n- 

In the following we study how the spin-dependent 
potentials influence the composite phase of the Mott- 
insulating core in the bulk. In Fig. [51 we show the 
GS density distributions of an unpolarized system of 
— = 20 under the influence of the different trap 
imbalances with Vj, — 6.0 x 10~^ and a strong repul- 
sive interaction of w = 6. For 7 = 1, a Mott phase is 
formed in the bulk region of the trap. With the increase 
of the confinement for the spin-up atoms, the spin-down 
atoms are repelled from the center of the trap. The 
Mott phase induced by the interaction between the lo- 
cally spin-balanced fermions is changed into the Mott- 
like phase induced by both the interaction between the 
locally spin-imbalanced fermions and the spin-dependent 
potentials. At the critical point of 7 = 2.6, the phase 
separation starts [see Fig. ^h)] and the strong confine- 
ment for the spin-up atoms forms an insulating core of 
fully polarized fermions in the center of the trap, over 
which the local occupancy is a unit. This insulating 
core is regarded as a ferromagnetic insulating phase since 
it is incompressible in nature and originates from the 
Pauli exclusion principle, which differs from the Mott- 
insulating phase induced by the repulsive interaction be- 
tween fermions [s^, such as in Fig. ^a). The unit 
core becomes stable for 7 > 6 by further increasing 7. 
Upon reaching the phase separation point and beyond, a 
plateau of constant rrii = 0.5 is formed in the center of 
the local magnetization. 

In Fig. [7] we show the case of a polarized system of 
= 30 and Ni = 10 with a strong repulsive interaction 
of M = 6. For 7 = 1, the spin- up fermions form Wigner- 
lattice-type profiles inside the Mott core which oc- 
curs at low fillings or equivalently at large u and can be 
explained by mapping Eq. ([l} into the antiferromagnetic 
Heisenberg model ^9|. We notice that, compared to the 
unpolarized system of that in Fig. [51 the polarized sys- 



5 





(a) y=1 



(b) T=1 .45 





FIG. 6: (Color online) Density distributions for spin-up and 
spin-down fermions together with their sum (the total GS den- 
sity) and difference (the local magnetization) plotted against 
the site with strong repulsive interaction {u = 6) in the 
spin-dependent external potentials. The system consists of 
Nf = 20 and iV^ = 20 fermions with = 6.0 x 10"=^. Three 
different ratios of confining potentials are shown: (a) 7 = 1, 
(b) 7 — 2.6, and (c) 7 = 6. 



tern becomes more easily reaches the phase separation 
point (7 = 1.45). That is, the Mott phase in the po- 
larized system is less robust against the increase of the 
interaction strength and confinement. 



IV. CONCLUSIONS 

In this article we perform a theoretical study of a ID 
Fermi-Hubbard model in a spin-dependent harmonic trap 
within the DMRG techniques. The interplay between the 
repulsive interaction and the spin-dependent harmonic 
trap is studied for the system of spin-balanced or spin- 
imbalanced Fermi gases. We find that, for the system 
in the spin-dependent external potentials, there exists a 



FIG. 7: (Color online) Same as Fig. [6] but for the spin- 
polarized system of — 30 and A^^ = 10. Three differ- 
ent ratios of confining potentials are shown: (a) 7 = 1, (b) 
7 — 1.45, and (c) 7 = 6. 



critical interaction strength beyond which a phase sep- 
aration can occur with two Fermi components staying 
in the different spatial regions. For the system with a 
weak interaction strength, upon increasing the trap im- 
balance, the spin- up atoms are confined more and more 
in the center of the trap and a depletion occurs for the 
spin-down atoms due to the intercomponent repulsive in- 
teractions. However, the weak repulsive interaction be- 
low a critical value is not capable of achieving a full phase 
separation. For the system with strong intercomponent 
repulsive interactions, a complete phase separation is re- 
alized at the strong confinement imbalance where spin- 
down atoms are repelled out of the bulk region with only 
spin-up atoms remaining. 

For the system with both strong confinement and 
strong repulsive interactions, where a composite phase of 
the Mott-insulating core is formed in the center, we show 
that, upon increasing the trap imbalance, the Mott phase 
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induced by the interaction between the locahy spin- 
balanced fermions is changed into the Mott-hke phase in- 
duced under the interplay between the interaction of the 
locally spin-imbalanced fermions and the spin-dependent 
confining potentials. Upon reaching the phase separation 
point and beyond, the ferromagnetic insulating phase due 
to the Pauli exclusion principle appears, which is of the 
unit core. In the distribution of the local magnetization, 
a step structure contributed by spin-up atoms alone is 
formed with a big slope from a negative to positive value. 
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